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Simulation  of  Matched  Field  Processing  for  Convergence  Zones 
in  a  Strong  Bottom  Field 
by  D.  F.  Gordon 

To  be  delivered  at  the  113th  Meeting  of  the  Acoustical  Society  of 
America,  Indianapolis,  Indiana,  Session  DD,  14  May  1987 

In  most  wave  theory  computations  of  underwater  sound  propagation  the 
convergence  zone  (cz)  only  field  can  be  computed  much  more  easily  and 
accurately  than  can  a  complete  field  that  includes  bottom  reflected  enr :  y. 

For  this  reason  the  effect  of  using  the  cz  only  field  on  matched  field 
processing  was  investigated.  The  investigation  was  extended  to  other 
simplified  descriptions  of  the  matching  field  including  simple  fields 
equivalent  to  alternating  polarity  arrays.  Some  results  of  this  investigation 
will  be  reported  here. 

At  frequencies  below  30  Hz  the  bottom  reflected  field  may  be  stronger 
than  the  cz  field  at  the  range  of  the  second  cz.  This  is  illustrated  in  the 
first  two  viewgraphs.  Viewgraph  1  shows  a  contoured  propagation  loss  field 
showing  the  first  and  second  cz.  Also  shown  are  the  depth  of  the  16  elements 
of  a  vertical  array.  The  spacing  between  elements  is  1/2  wavelength  at  30  Hz. 

We  will  be  interested  later  on  in  the  pressure  amplitude  along  the  array. 

Note  the  indication  here  that  there  will  be  alternating  intervals  of  high  and 
low  amplitude  down  the  array. 

The  contoured  field  we  see  here  is  for  a  100m  deep  source  and  the  array 
is  at  100  km  range.  Sound  pressures  shown  here  are  computed  by  a  normal  mode 
program  for  a  Pacific  profile  at  30  Hz  frequency.  A  very  lossy  bottom, 
modeled  as  a  negative  gradient  half  space,  was  used  to  suppress  bottom 
reflected  energy. 

Viewgraph  2  shows  a  similar  computation,  but  here  a  realistic  bottom  was 
used.  The  bottom  was  modeled  as  two  fluid  layers  representing  70m  of  fine 
grained  sediment  over  rock.  Note  that  the  convergence  zones  can  still  be 
distinguished,  but  are  distorted  by  the  bottom  reflected  energy. 

The  array  at  100  km  range  is  again  shown.  The  16  pressures,  expressed  as 
complex  numbers,  for  this  source  depth  of  100  m  and  this  range  will  be  used  as 
the  target  vector.  The  16  pressures  for  each  point  on  the  total  grid  of 
source  depths  and  ranges  constitutes  the  field  vectors. 

Let  us  first  look  at  the  results  of  simulating  matched  field  processing 
for  this  total  field  case.  The  total  grid  comprises  the  range  points  and  10 
different  source  depths  from  20  to  200m.  Viewgraph  3  shows  the  result  of  this 
simulation.  The  contours  are  the  squared  absolute  value  of  the  correlation 
coefficient  between  the  16  pressures  at  the  target  vector  and  any  field 
vector.  The  expected  high  correlation  can  be  seen  at  100m  source  depth  and 
100m  range.  However,  there  are  many  secondary  peaks  nearby,  which  in  the 
presence  of  noise  might  be  mistaken  for  the  primary  peak. 

We  will  now  inquire  into  the  effect  of  using  a  cz  only  field.  This  will 
be  done  by  substituting  the  sound  pressures  computed  in  this  way  for  the 


target  vector.  The  best  simulation  of  reality  would  require  that  the  numerous 
field  vectors  rather  them  the  target  vector  be  re-computed  by  the  method  being 
considered.  However,  to  minimize  computing,  the  field  vectors  are  always  the 
same,  computed  from  the  total  field.  It  is  assumed  that  results  for  the  two 
methods  will  differ  in  detail  but  be  similar  in  general  behavior. 

Viewgraph  4  repeats  viewgraph  3  but  for  an  enlarged  view.  Here  we  see 
the  true  target  point  at  100  km  range  and  also  a  second  area  of  enhanced 
correlation  at  10  km  greater  range.  This  turns  out  to  be  a  location  where  the 
upgoing  rays  at  the  target  have  refracted  and  are  returning  downward  in  a 
similar  pattern.  The  cz  only  target  vector  correlation  will  now  be  shown  on 
this  same  scale. 

Viewgraph  5  shows  this  matched  field.  Note  first  that  we  no  longer  have 
a  100  percent  correlation  point.  In  fact,  the  target  location  no  longer  is 
the  strongest  correlation.  This  occurs  for  about  80m  greater  source  depth. 
This  correlation  function  is  much  smoother  and  more  massive  than  the  previous 
full  field  one.  However,  it  no  longer  has  the  source  localization  properties 
of  the  previous  one.  Again  we  see  the  downgoing  part  of  the  cz  about  10  km 
farther  in  range. 

Let  us  now  breifly  look  at  the  correlation  process.  Viewgraph  6  shows 
the  correlation  equation.  Here  t  is  the  target  vector  and  p  is  the  field 
vector.  The  summations  are  over  the  16  array  elements.  Note  that  t  can  also 
be  thought  of  as  the  weights  and  delays  for  an  ordinary  plane  wave  beam 
former.  The  value  of  the  correlation  is  then  the  response  of  the  beamformer 
to  a  source  at  the  source  depth  and  range  of  the  field  vector.  In  the 
remainder  of  this  paper  the  target  vector  will  be  viewed  interchangeable  in 
these  two  ways. 

Viewgraph  7  shows  the  sound  pressure  from  element  to  element  down  the 
array.  The  complex  pressure  is  shown  as  an  amplitude  and  phase.  The  left 
hand  plot  is  the  pressures  of  the  cz  only  target  vector  used  previously.  Note 
the  alternating  high  and  low  amplitude  down  the  array.  There  is  a  tendancy 
for  the  phase  to  advance  180°  from  one  amplitude  peak  to  the  next.  The  right 
hand  plot  is  the  pressure  that  would  be  set  up  by  two  equal  plane  waves 
intersecting  the  array  from  10°  above  and  below  the  horizontal.  The 
interference  pattern  set  up  down  the  array  is  a  sine  wave.  A  comparison  of 
these  two  plots  indicates  that  the  cz  pressures  have  a  strong  component  of  two 
plane  waves  near  +  10° .  The  overall  amplitude  of  the  target  vector  is 
normalized  out  of  the  correlation  function  so  it  is  only  the  relative 
amplitude  and  phase  that  matter  in  the  target  vector.  However,  the  location 
of  the  sine  wave  pattern  in  depth  along  the  array  is  a  significant  parameter 
and  depends  upon  the  relative  phase  of  the  two  interfering  plane  waves.  Here, 
their  relative  phase  has  been  chosen  to  match  the  computed  target  vector. 

Viewgraph  8  shows  two  more  simplified  target  vectors.  The  left  plot 
represents  steering  to  a  single  plane  wave  arriving  from  10°  below  the 
horizontal.  The  right  plot  represents  an  alternating  polarity  array.  This 
steering  can  be  achieved  by  connecting  sequences  of  array  elements  in 
alternating  polarity.  Such  steering  is  of  interest  in  matched  field 
processing  because  the  multiplication  in  the  correlation  process  can  be 
omitted  and  field  vector  points  are  only  added  or  subtracted.  This  greatly 
reduces  the  computation  load.  The  alternations  of  polarity  here  are  chosen  to 


match  those  of  the  previous  viewgraph. 

The  correlation  function  over  the  field  will  now  be  shown  for  these 
target  vectors.  Viewgraph  9  shows  contours  for  matching  to  a  10°  upcoming 
wave.  The  0.4  contour  here  covers  about  the  same  area  as  the  same  contour  did 
in  the  cz  matched  field  shown  earlier.  The  0.2  contour  here  spreads  out 
substantiably  farther  than  did  the  similar  one  of  the  cz  case.  Note  that  very 
little  of  the  downgoing  part  of  the  cz  correlates  with  this  target  vector. 

Viewgraph  10  shows  the  0.2  contour  only  for  the  target  vector  created  by 
±  10°  plane  waves.  This  contour  is  much  narrower  than  the  previous  for  a 
single  10°  plane  wave,  but  the  two  sections  of  the  cz  are  almost  equally 
correlated. 

Finally,  viewgraph  10  for  the  alternating  polarity  steering  shows 
correlations  that  are  similar  but  not  as  sharp  as  the  previous  case.  An 
application  of  this  concept  would  be  observing  only  the  relative  phase  down  an 
array,  quantizing  to  +  and  -  and  then  adding  or  subtracting  the  computed  field 
vectors  to  obtain  the  correlation. 

In  conclusion:  Simulations  show  that  if  field  matching  at  30  Hz  were 
attempted  with  a  field  computed  only  from  cz  energy,  correlation  functions 
would  be  broad  and  smooth,  but  would  not  peak  at  the  true  source  depth. 
Further,  secondary  peaks  would  be  found  at  locations  in  the  cz  where  arrival 
patterns  are  in  some  way  similar. 

It  is  shown  that  results  quite  similar  to  those  obtained  by  using  the  cz 
field  could  be  obtained  by  methods  similar  to  plane  wave  steering  and  even 
alternating  polarity  steering. 
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